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Abstract 

We study social cost losses in Facility Location games, where n selfish agents install facilities 
over a network and connect to them, so as to forward their local demand (expressed by a non- 
negative weight per agent). Agents using the same facility share fairly its installation cost, but 
f-H ■ every agent pays individually a (weighted) connection cost to the chosen location. We study 

I the Price of Stability (PoS) of pure Nash equilibria and the Price of Anarchy of strong equilibria 

(SPoA), that generalize pure equilibria by being resilient to coalitional deviations. A special case 
Q I of recently studied network design games. Facility Location merits separate study as a classic 

model with numerous applications and individual characteristics: our analysis for unweighted 
agents on metric networks reveals constant upper and lower bounds for the PoS, while an O(lnn) 
^ ' upper bound implied by previous work is tight for non-metric networks. Strong equilibria do 

CNJ I not always exist, even for the unweighted metric case. We show that e- approximate strong 

0^ ■ equilibria exist (e — 2.718. . .). The SPoA is generally upper bounded by O(lnW) (W is the 

sum of agents' weights), which becomes tight O(lnn) for unweighted agents. For the unweighted 
metric case we prove a constant upper bound. We point out several challenging open questions 

' that arise. 

O " 
00 

'. 1 Introduction 
> . 

^ . Modern computer and communications networks constitute an arena of economic interactions 

^ I among multiple autonomous self-interested entities (network access providers, end-users, electronic 

commerce enterprises, content storage and distribution enterprises). The internet is perhaps the 
most massive and global field where economic networking interactions take place. The recently 
established study of network formation games [14] (chapter 19) aims at understanding how the 
competitive (or coalitional) activity of multiple such selfish agents affects the network's character- 
istics and its efficiency. In this paper we consider the setting of n Content Distribution Network 
enterprises interacting over a network. Enterprise i takes decisions on where to store replicas of 
digital content over the network, so as to satisfy access demand of local customers, situated at node 
Ui. Demand is expressed by a non-negative weight wt for enterprise i. Every enterprise chooses 
strategically a location v on the network for instahation of content replicas, so as to minimize 
its individual expenses for (a) storage/management of the content, and (b) weighted connection 
(bandwidth/delay) costs to the chosen location. We study a Facility Location game played among 
selfish agents (enterprises), that naturally models this situation. We use a fair cost allocation rule 
- known as Shapley cost-sharing [3] - for facility installation costs: for a location v installation cost 
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Pv is shared in a fair manner among agents that receive service from v, so that agent i pays , 
where W(v] denotes the total demand forwarded to v. The weighted connection cost is given by 
Wid(ui,v), where d(ui,v) denotes connection cost per unit of demand. 

Our work focuses on bounding the social cost of cost-efHcient and stable network infrastructures; 
stable networks will be represented by pure strategy Nash equilibria and strong equilibria of the 
corresponding Facility Location game. The social cost will be the sum of individual costs experi- 
enced by the agents-players of the game. Strong equilibria are an extension of pure Nash equilibria, 
dicussed in [2, 9], and essentially introduced by Aumann in [5]. Strong equilibria are resilient to 
coalitional pure deviations: a strategy profile is a strong equilibrium if no subset of agents can 
deviate (by jointly adopting a different pure strategy), so that all of its members are better off. 
Thus this notion captures the possibility of coalitional behavior among agents. Occurence of such 
behavior is naturally expected in rapidly evolving modern markets, especially the ones involving 
the exploitation and exchange of digital goods and services. 

We derive bounds on the Price of Stability (PoS) of pure equilibria, defined as the cost of 
the least expensive equilibrium relative to the socially optimum cost [3]. For strong equilibria we 
derive bounds on their (strong) Price of Anarchy (SPoA), i.e. the cost of the most expensive strong 
equilibrium relative to the socially optimum cost. Let us note that the Price of Anarchy of pure 
Nash equilibria was introduced in [12] as the cost of the most expensive pure Nash equilibrium 
relative to the socially optimum cost, and measures essentially social cost losses incurred due to 
selfishness of agents, and lack of coordination among them. The PoS on the other hand measures 
social cost losses only due to selfishness: if all players coordinate (even by following some externally 
provided instructions, or by interacting on the basis of a common protocol) they may reach the 
least expensive equilibrium. The notion of strong equilibria inherently permits coordination among 
subsets of agents. Therefore we can intuitively expect the SPoA to almost match the PoS. Our 
results confirm this intuition. 

The study of the price of stability in network design games with fair allocation of network 
link costs was initiated in [3]. One challenging problem remaining open since then is improving 
upon an upper bound of PoS = O(logn) and a lower bound (recently shown in [10]) of -y, for 
unweighted players on undirected networks. A series of recent works [10, 6, 7, 9, 1] provides results 
(polylogarithmic upper bounds) with respect to the social cost of pure Nash and (approximate) 
strong equilibria in the network design game model of [3], also for the case of weighted players. We 
review these results in section 2. We also explain how the Facility Location game that we study 
is a special case of the more general model introduced in [3], that also includes "delay" costs in 
using network links (refered to as connection costs in the case of Facility Location), apart from 
fairly allocated "installation" costs. For the case of metric connection costs we were able to prove 
greatly improved and almost tight bounds for the PoS and the SPoA; this makes the metric Facility 
Location game and exceptional special case of the model introduced in [3] . 

Summary of Results We analyze the PoS of the unweighted metric Facility Location game 
(section 4), and prove constant upper and lower bounds. Note that an O(lnn) general upper 
bound implied by the work of [3] is tight for non-metric networks (we discuss this in section 3). 
Our technique relies on direct analysis of social cost evolution during an iterative best response 
performed by the agents, until they reach equilibrium. Wc show in particular that given any 
initial configuration (including the social optimum), the social cost of the reached equilibrium 
is at most 2.36 times the initial social cost and at least 1.45 times this cost in the worst case. 
Metric Facility Location is the first case of the model proposed in [3], in which an additional 
structural network property (triangle inequality) yields a constant almost tight Price of Stability. 



2 



Metric 
Unweighted 



Unweighted 



Non-metric 
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PNE I SE 
PoS 
SPoA 



V I e-apx 
G (1.45,2.36) 
0(1) 



a/ I e-apx 
G(lnn) 
G(lnn) 



e-apx I e-apx 



O(lnW) 
O(lnW) 



Table 1: Summary: Facihty Location Games with Fair Allocation of Facility Costs. 

It is also interesting that incorporation of link delays (connection costs) makes the previously 
derived 0(H(n)) upper bounds tight for the non-metric unweighted game. Subsequently we study 
strong equilibria (section 5). Although they do not always exist, we prove that a-approximate 
strong equilibria exist (no subset deviation causes factor a improvement to all of its members), 
for a > e = 2.718... in general networks and weighted agents. For the Price of Anarchy of 
approximate strong equilibria (SPoA) we show an O(lnW) general upper bound which becomes 
Gflnn) for unweighted agents (section 6). For the metric unweighted case we prove a constant 
upper bound (section 6.1). 

Except for the SPoA analysis in paragraph 6.1 (theorem 4), the rest of the described results 
have appeared in [11], in a shorter version ^. We note that the analysis of the SPoA of approximate 
strong equilibria for the metric unweighted case in paragraph 6.1 essentially generalizes the analysis 
of the corresponding result appearing in [11] (theorem 2); this result concerned the SPoA of exact 
strong equilibria only when they exist. 



Anshelevich et al. first studied the Price of Stability for network design games with general [4] and 
fair cost allocation [4]. In the latter case n agents wish to connect node subsets over a network, 
by strategically selecting links to use. Every network link is associated to two cost components, an 
installation cost and a delay cost. Link installation cost is shared fairly among agents using the 
same link in [3], while every agent experiences a delay cost given by a polynomial function of the 
number of agents using the link. The authors showed that for unweighted agents these network 
design games belong to the class of potential games introduced by Monderer and Shapley in [13], 
hence they have pure strategy Nash equilibria. In particular, a potential function defined in [13] is 
associated to these games in the following manner: an arbitrarily initialized iterative best response 
procedure carried out by the players reaches a local minimum of the potential function, which 
corresponds to a pure strategy Nash equilibrium of the game. In [3] the authors developed elegant 
arguments using the potential function, so as to upper bound the PoS for unweighted agents. In 
case of polynomial delay costs of degree at most k the PoS was shown to be at most 0((k-|- 1 ) Inn). 
For the case of directed networks this bound was shown to be tight (for k=0 and zero delays), 
whereas for undirected networks a lower bound of | was given. This was recently improved to 
^ by Fiat et al. [10], who also showed an upper bound of O (log log n) for PoS, when exactly one 
agent resides on each node of the network and all agents need to connect to a common sink node 
(single-sink case). For weighted agents a potential function proving existence of pure equilibria was 
developed for only 2 agents in [3] . Chen et al. [7] showed that pure strategy Nash equilibria do not 
always exist for weighted network design games with fair allocation of link installation costs. They 

result of 2.36-approximate strong equilibria that appeared in a previous version of this report was erroneous 
and has been removed. 
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developed a tradc-off for the social cost of approximate equilibria versus the approximation factor, 
as a function of the maximum weight taken over all agents. 

An extensive study of strong equilibria for network design games with fair cost allocation (but 
without delays) appeared recently by Albers in [1], for unweighted and weighted agents. Strong 
equilibria do not always exist, but she showed that 0(H(n)) and 0(ln W]-approximate strong 
equilibria do exist (W is the sum of the agents' weights) for unweighted and weighted games 
respectively. She proved 0{H(n)) and O(lnW) upper bounds for the SPoA of approximate strong 
equilibria which are tight for directed graphs. For undirected graphs she showed Cl{^/logn) and 
£l(\/logW) lower bounds. Finally and most importantly, she showed an ^( iogk)gw ) lower bound 
for the PoS of weighted network design games. Strong equilibria in the context of (single-sink) 
unweighted network design games with fair cost sharing were first studied in [9] . The authors gave 
topological characterizations for the existence of strong equilibria on directed networks, and proved 
that SPoA = e(logrL). 

As of the recent literature, the gap for the identification of the PoS of network design games 
with fair cost allocation remains open for unweighted agents, with the lower bound of y= being 
the best known so far, to the best of our knowledge. In effect, this means that the impact of an 
undirected network structure to the PoS is not well understood, since the upper bounding potential 
function arguments developed in [3] do not incorporate network structure. The Facility Location 
game is a special case of the model studied in [3], that is interesting on its own right: it finds 
numerous applications and exhibits intriguing characteristics. It emboddies non-shareable delays 
explicitly and in a sense specializes single-sink network design considered in [9, 7]: we can simply 
augment the network with a node t and set links (v,t) to have fairly shareable cost, equal to the 
facility opening cost at v. The original network links have a delay cost only. Then every agent 
needs to choose at most two edges from the node it resides on, to t. 

Facility Location Games An unweighted metric facility location game with uniform facility 

costs was studied in the context of selfish caching in [8]. The network model of that work is 
essentially equivalent to the one we discuss here, apart from the fact that fair cost-sharing of 
facility costs was not used: an agent could connect to a facility payed exclusively by another agent. 
Another difference is that facility opening cost was uniform accross all nodes of the network, whereas 
we consider node-dependent costs. The PoA and PoS of this model were shown to be unbounded. 
The authors devised an extension of their game model, with payments exchanged among agents, 
in which the socially optimum configuration of the original game is rendered a pure strategy Nash 
equilibrium (hence PoS = 1 for the extended game). Vctta studied a class of games for competitive 
facility location [15], for which he proved existence of pure strategy Nash equilibria and an upper 
bound of 2 for the PoA. In this competitive setting enterprises open facilities at certain nodes 
and try to attract customers to connect to them. He also illustrated applicability of his model 
in the context of the k-median problem model. In [6] the authors considered single-sink network 
design with fair cost-sharing of all resources (network links) used by agents. They argued how this 
model can be viewed as a Facility Location model, but did not take connection costs (non-shareable 
delays) into account. 

3 Definitions and Preliminaries 

The network will be a complete graph G(V, E), each edge (u, v) of which is associated to a non- 
negative cost d(u,v). We consider a set A of n agents. Each agent i resides on a node iXi G V, and 
is associated to a non-negative demand weight Wi. The strategy space of agent i is V: i chooses 



4 



a location v G V, where it will receive service from. Opening a facility at node v G V costs (3v. A 
strategy profile (configuration) is denoted by s = (si, . . . .Sn), where Si G V. The total weight of 
agents receiving service from v under s is Ws(v). The cost Ci(s) experienced by agent i in s is: 

Ci(s) =Wid(iii,Si) + 



Every agent i pays a fraction y^^^. ^ of the facility installation cost at st. Denote by Fj C V the set 
of facility locations specified under s. The social cost c(s) is then: 

C(s) = Y_ = X "^i^^^i- Si) + JL WJs^ ^ Wid(Ui, Si) + ^ |3v 

We use W(I) for the sum of weights of agents in set 1. By Ci(s) we denote the total cost of agents 
in I under s. Furthermore, given a facility v G Fs) we use Cv(s) to denote ^i-g. =v ''i-(s)- 

Pure Nash Equilibria The unweighted Facility Location game specializes network design games 
first studied in [3], and is a potential game [13]; it is associated to a potential function, that can be 
locally minimized by iterative best response performed by players, and the reached local minimum 
corresponds to a pure Nash equilibrium. In iterative best response we choose iteratively an arbitrary 
player i, and let him/her decide a strategy that minimizes its individual cost with respect to the 
current configuration s^i = (si , . . . , Si_i , Si+i , . . . , s^) for the rest of the players. The process ends 
when no player i can improve his/her individual cost under s_i. Then s is a pure Nash equilibrium. 

A simple example shows that the PoA of pure equilibria can be n in the worst case. Take a 
network of two nodes u, v. Opening a facility on either node has a cost |3. We set d(u, v) = 
and let n agents reside on u. Then if every agent i plays St = v, i = 1 . . . n, the resulting 
configuration is a pure Nash equilibrium of cost np, because every agent pays exactly |3 and has 
no incentive to open a facility at node u by paying the same cost (3. The social optimum consists 
of a single facility opened on u and all agents being locally serviced by it, at zero connection cost. 
Thus the socially optimum cost is |3. 

Potential function arguments developed in [3] yield an upper bound of H(rL) for the PoS, where 
F-l('n.) is the n-th harmonic number. This bound is tight for unweighted non-metric Facility Location. 
Fig. 1 presents a non-metric lower bounding example for unweighted agents. Assume uniform 
facility costs equal to 1. 2n unweighted agents reside on boldly drawn nodes, n of them on vj- 
The social optimum consists of facilities on vi and Vz, where Vi serves only the n agents residing 
on V2- The rest are served by V) . However they deviate to V2 from vi one by one. It is then: 
PoS > (2H(n) — H{2n)) / [2 + ne) > (Inn — ln2 — l)/(2-|-ne) = O(lnn). When it comes to metric 
networks however, the PoS is constant as we show in the following section. This result essentially 
indicates the impact of the network's structure (triangle inequality) to the social cost of efficient 
pure Nash equilibria. 

Definition 1 (Strong Equilibria) [2, 5] A strategy profile s is a strong equilibrium if no subset 

of agents can deviate in coordination, so that each and every one of its members is better off. It 
is an oi- approximate strong equilibrium, if no subset of agents ca,n deviate in coordination, so that 
each and every one of its members is better off by a factor strictly more than a > 1 . 



4 Unweighted Metric Facility Location: The Price of Stabihty 

We analyze evolution of an equilibrium through iterative best response, initialized at configuration 
s*. In fact, for any strategy profile s that is reached through iterative best response initialized at s*. 
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Figure 1: Lower bounding example for unweighted (non-metric) PoS. 



we will upper bound the ratio c(s)/c(s*). Thus we do not require that s* be the socially optimum 
configuration, nor that s is an equilibrium strategy profile. We only require that s is reached by 
iterative best response initialized at s*. Then taking s* to be the socially optimum configuration 
and s to be the reached equilibrium will yield an upper bound for the Price of Stability, as a 
corollary. To facilitate clarity, we give a brief description of the analysis' plan. Given any facility 
node V G Tg* denote by As*(v) = {i|s^ = v} the subset of agents connected to v in s*. Then Cv(s*) 
will be the total cost incurred by As*(v) collectively. We will upper bound c(s)/c(s*) as follows: 



c(s*) " LvGF^* LieA3.(v)Ci{s*) ~ LveF3«Cv(s*) -™ Cv(s*) 

Define Ac(i) = Ci(s) — Ci(s*) to be the increase in cost caused by agent i during iterative best 
response initialized at s* and reaching s. For any subset A' C A of agents we also use Ac(A') = 
^^g^/(Ci(s) — Ci(s*)). To evaluate the upper bounding expression (1), we are going to use an upper 
bound on Ac(As*(v)) = Hvs?^v^'^(^)' valid for any v G Fj*. This is the increase in social cost 
caused during iterative best response by agents connected to v in s*. Then for any v G Yg* we 
will maximize '^"(s by lower bounding Cv(s*) for any v G Fg*. In determining an upper 

bound on Ac(As*(v)) we prove the following lemma, which charges any specific agent i a bounded 
amount of social cost increase during the execution of iterative best response. 

Lemma 1 Let Ag* (v) be the subset of agents that are connected to v in s* . For any i G Ag* (v) 
define A^* (v) C Ag. (v) to be the subset of agents that have not yet deviated from v, exactly before 
the first deviation of i. Then Ac[i) < (Sv/lAg, (v)|. 

Proof. For simplicity let lA^. {v)| = lCi(v). Clearly i G A^, (v). Let us analyze contribution of i to 
social cost increase during its first deviation. By deviating, i reduces its individual cost from Ci to 
c|, by joining another facility node v'. Then: 

Ci = Xi(v) + -— ^, C^ = Xi(v) + 



UvV ' ' At(v') 

Xi(v) and Xi(v') is the connection cost payed by i before and after its first deviation. Ai(v') is the 
number of agents sharing facility cost at v', including i. Since Ci(v') < Ci(v), it is Xi(v') — Xi(v) < 
icfe ~ \(v') ' ^Cv(i) = c[ — Ci be the difference caused to the social cost by this single deviation 
of i. We examine the following cases: 
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1. ki(v) > 1, Ai(v') > 1: Then ACv(i) = Xi(v') - Xi(v) < - xf^- 



2. ki(v) = 1 , Ai(v') > 1 : Then Acv(i) = -|3v + Xi{v') - Xi(v) < 

3. ki(v) > 1, Ai(v') = 1: Then ACv(i) = Pv' + Xi(v') - Xi(v) < 

4. ki(v) = 1, Ai(v') = 1: Then ACv(i) = |3v' - |3v + Xi(v') - Xi(v) < 0. 

Clearly the above hold in general during execution of iterative best response, for any agent that 
performs a single deviation from a node v to a node v'. 

Now we implement a charging procedure along with iterative best response. Give all agents an 
initial label = i, before executing iterative best response. The labels will be updated during 
iterative best response, so that increases caused by agent i will be charged to an appropriate agent 
Initialize Ac(£(i)) = for all i G A. We will make use of the auxiliary variables ki(v) and Ai(v'), 
as they were described before. At any time Ai(v) is the number of agents (including i) connected to 
V, right after i joined v. We need to initialize At(v) to a value for every v G fg* and i G Ag* (v), to 
implement the charging scheme. Set = A| to a distinct value from {1,2,..., lAs*(v)|}. ^ni)iv) 
will be always (during iterative best response) the number of agents connected to v exactly before 
deviation of i from v (before i leaves v). Charging is then implemented by relabeling deviating 
agents in the following manner. 

1. If k£{t)(v) = A£(|)(v) no relabeling is needed. 

2. Otherwise there must be some j 7^ i connected to v such that A£(j)(v) = k^^i^iv). In this case 
exchange labels of i and j. 

Subsequently add the increase caused by deviation of i to Ac(€(i)). Finally, set X^[i){v') equal to 
the number of agents connected to v' right after i has joined v'. 

By the previous definitions it follows that if k^^i) (v) / X^^^ (v) , then it is always k^f^j (v) > A^fi^) (v) , 
i.e. i has joined v before some agent j with A£(j)(v) = k£(i)(v), but leaves v "out of order", i.e. before 
j leaves. By exchanging labels of i,j we add the increase caused by i to the agent that previously 
labeled j. Possible increases in l.,2.,3.,4., imply that any agent is charged by the end of iterative 
best response at most ^^t'^'j^^^ for some i. If we "guess" the exact order of first deviation of all 

agents, we can initialize Ai(v) = ki(v). Then each i itself will be charged at most 1..^", □ 

y, Ci(s) 

In what follows we are going to upper bound the ratio — ^^^"(s*') — ^^i any v G Fg* . If no agent 

II- ci(s) 

i G As* (v) ever deviates from playing v, then clearly it will be — ^'^'^"(s*'] — ~ ^ • analysis 
focuses on two remaining cases: either (i) there is a non-empty set of agents As(v) C As*(v) that 
never deviate from v during iterative best response or (ii) all agents i G As* (v) deviate from v during 
iterative best response. Case (i) is examined in proposition 1, whereas the analysis is concluded by 
analysis of (ii) in theorem 1. 

Proposition 1 Let s be a strategy profile reached by iterative best response initialized at strat- 
egy profile s*. For any facility v G Fg* define As*(v) = {i|s| = v} and let As(v) C As*(v) the 
subset of agents that never deviated from v, during iterative best response. If As(v) 7^ 0, then 
LieA^.M Ci(s) <2.36-Cv(s*). 
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Proof. As discussed previously, we will upper bound the ratio '^^^^ ^c^^^s*'^^*^^^^ ' deriving an 
upper bound for Ac(As*(v)) and a lower bound for Cv(s*). By lemma 1 follows immediately: 

Ac(As.(v)) < |3v[H(|A,.(v))-H(|AJv)l)] (2) 

because agents in As{v) never deviated, whereas the rest caused a cost increase equal to pv times 
a harmonic series term each. In what follows we show a lower bound on Cv(s*). 

Fix any facility v G Fg*, and define Is*(v) = As*(v) \ As(v). Without loss of generality assume 
such an order of agents, that agents i € Ig* (v) (and for every v G Fg*) "best-respond" consecutively; 
e.g. assume that the algorithm scans the facilities in Fj* in an arbitrary fixed order and for each 
facility, the agents connected to it in an arbitrary fixed order. We focus only on the first deviation 
of each agent i G Ig* (v) for some v G Fj*. Let v' be the node that i deviates to. For convenience we 
use xT = d{ui,v) and let 6x? = d(u|,v') — d(ui,v). Let be the total number of agents connected 
to v' right after deviation of i. The new cost of i right after its first deviation is: x| + 6x| + 
For any other agent j G Aj. (v) that, either deviates from v to some node v" after the deviation 
of i, or never deviated from v (in this case consider a trivial deviation with v" = v), we have 
respectively: 

d(uj,v)+6x; + %^<d(u^,v') + ^ (3) 
Aj Ai 

Substitute d(uj,v') in (3) by triangle inequality: d(uj,v') < d(uj,v) + d(ui,v) + d(ui,v'). Let 
k| denote the number of agents connected to v right before deviation of i from v. Also note that 
6x| + < because i decreases its individual cost by deviating. Thus we obtain: 

Aj A| / - 2 V ' ^ Aj k? 



d(ut,v) > 5x^ -S.l+-f-^)> + ^-^) (4) 



Because d(ut, v) > 0, and because the latter has to hold for every pair of distinct agents i, j G As* (v) 
we deduce: 



d(ui,v) > ^ 




|3v 
k? 



(5) 



We use (5) for the connection cost of agents in Is*(v), and for agents in As(v) = As*(v) \ 
Is*(v); in essence for every i G As(v) we simply set in (5) v" = v, k| = |As(v)|, and take 
maxj:s*=v (Sx? + = j^^f^i' The cost Cv(s*) is then: 

„ IAsN.(v]l 

c,(s-) > 3.+ ^ ^ (^-1) 

k=[As(v)| ' ^' 

r, , Pv/|As*(v)| - |As(v)| \ 

= Pv+y(^ \A~(vj\ H(|As*(v)|) + H(|As(v)|jj (6) 

Using the bounds (6) and (2) we deduce the following upper bound: 

\ I + IS'""" + H(|A..(v)|) -HdAsWP) 

.t:<v)"'7 ' ,+ .( IA..(;^HA.HI _H(|A..(v)|) + H(|A.(v)|)) 
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^ 1 .5 + |As* (v)|/|As(v)l + ln(|As» (v)|/|As(v)|) 
- 0.5 + |As* (v)|/|As(v)| - ln(|As* (v)|/|As(v)|) 



(7) 



The latter inequality follows by bounds for the harmonic number, i.e. y + lnm < H(m) < 1 +lnra, 
where y > 0.5 is the ii^u/er constant. By setting y = |As* (v)|/|As(v)|, we can numerically maximize 
the upper bound with respect to y, to at most 2.36, for y ~ 2.47. □ 

Theorem 1 If s is a strategy profile for the unweighted Metric Facility Location game, reached by 
iterative best response initialized at a strategy profile s* , then c(s*) < 2.36c(s). 

Proof. We only need to complement the result of proposition 1 by considering the case where 
As{v) = 0. That is, Is*{v) = As*(v), and all agents connected to v in s* deviate during iterative 
best response. If for each i G Is*(v) that deviates from v to v' (v' is not necessarily the same for 
every i) it is 6xt + < then: 

Y_ Ci(s)= Y. + + <Pv+ Y <=Cv(s*) 

Thus, assume there is at least one agent j with bx^ + ^ > -jXTM- general, let r be the largest 
integer, so that 5x? + ^ < ^ for all i G Is*(v) = As*(v). Using the same arguments that led 
to (5), and taking maxi + = we deduce: 

xt = d(Ui,v)>l(^^--^j, fork?>r and x? = d(iH,v) > 0, for k? < r (8) 

Using (8) we end up with a similar lower bound to (6) for 1 < r < |As*(v)|: 

lc=r 

To finish the proof, we note that by (8), deviation of agents with k? < r, yields a total cost increase 
(payed as new connection cost) of at most r x = p^, which is exactly equal to the cost saved by 
closing the facility at v (because all agents deviate). Thus we can simulate the situation witli the 
case analyzed in proposition 1, by setting As(v) = {i G As*(v)|k| < r} and r = |As(v)|, so as to end 
up with (7). This is technically equivalent to assuming that agents i with k| < r never deviated. □ 

Corrolary 1 The Price of Stability for the unweighted Metric Facility Location game with fairly 
allocated facility costs is at most 2.36. 

4.1 A Lower Bound on the Price of Stability 

We analyze a worst-case example, that makes the PoS for unweighted metric Facility Location at 

least 1 .45 asymptotically. Experimental evidence shows a lower bound > 1 .77, which we believe is 
tight for PoS. This is because the analysis of theorem 1 embodies some losses due to bounding of 
harmonic numbers by logarithms. 

Our construction appears in Fig. 2(a). Take 2n agents, n of them residing on a singe node v of 
the network. Facility opening costs are assumed 1 everywhere. The social optimum s* has 1 + ^/r^ 
facilities: v and v]^, I = 1 . . . k, where k = -^/n. Let the n agents residing on v be serviced by v in 
s* , whereas the rest are equally partitioned to facilities V| , I = 1 . . . k; every facility services a 
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n 



PoS 




(a) Lower bounding example. 



10^ 

1 \J 


1 ^2471 




1 69106 


1 \J 




10^ 


1 .76367 


10^ 


1 .76927 


107 


].77m 


108 


M7\6 



(b) Obtained estimates. 

Figure 2: Experimental estimation of a PoS lower bound for the unweighted metric case. 



batch of k = yri agents. In the least expensive equilibrium s, all agents from each facility of s* 
will deviate and miss-connect to v. We will analyze for any single facility - henceforth denoted 
by V* - the ratio (^tg^ *(v*) ^^iC^)) /Cv*(s*). Analysis for the rest of the facilities will be similar, 
and will yield asymptotically the same result. In particular, by choosing to analyze any batch of k 
agents that deviate from v* to v, the PoS will be lower bounded as: 



PoS 



lim 

n— >oo 



Ci(s] 



> lim 



ieA3*(v*) 



Ci(s) 



n-^oo (1/Vu + Cv«(s*)) 



(10) 



For some constant p G (0, 1 ) we will have r = |"(1 — p)k] of the agents in As* (v) increase their 
connection cost significantly by deviating, hence the social cost. For every i € As*(v) let be the 
distance of i from v*, and Sx? the increase in connection cost after deviation. We determine x? and 
6x| by following best responses of agents in As*(v*). Let A > n be the number of agents already 
connected to v before agents of As*(v*) deviate to v in order. In particular: 

• For i = 1 . . . r set 6X| = — ^ — e; a small decrease e > causes i to deviate to v. 

• For each of the remaining k — r agents set 6x* = ^^^rp^ ■ 

By similar arguments as for the derivation of (3), (4) and triangle inequality, we obtain x? = 
d(iii,v*) = max{0, ^ImaxSx? — 6x|)}, where max^ Sx? = This yields x| = for i > r (the 

remaining k — r agents). Note that the k — r remaining agents will deviate to v, because they prefer 
to pay strictly less than + instead of a cost share each, for the facility at v*. By 

deviation they "shut-down" the facility at v* and decrease the social cost by 1 , but pay a total 
connection cost ^^^^ for v, which tends to 1 for large n. Let us now derive the cost Cv*(s*). For 
convenience we use the notation AH(p, q) to denote H(p) — H(q). For the first r agents, summing 
up as in (6) yields: 



1 + 



1 



k-r + 1 



AH(k,k-r) 



77^ + 7AH(A + r,A) 
2 A + r 2 



The cost of agents in As*(v*) after deviation is: ^igA *{v) Ci(s) = Cv*{s* 
By substituting ^[^^ 6xT = AH(k, k-r) - AH(A + r, A), we obtain: 



' ^ k-r+l ^ 2-i 



(11) 



=1 
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iGAs* (v) 



k-r 1^ 
k-r+1 ^2 



k-r+1 



+ AH(k,k-r) 



2 VA + r 



+ AH(A + r,A) 



(12) 



We simplify (11) and (12) and substitute to (10) appropriately. For AH(k, k— r) we use the following 
bounds: 

y-1+ln 1— < AH(k,k-r) < 1 -Y + ln ^ 



where y > 0.5 is the Euler constant. Also we use ^ — k— r+l — "^p^' because (1 — p)k — 1 < 

r= [(1 -p)k] < (1 - p)k+1. Then (10) becomes: 



p-lA 



1-p 



PoS > liin 



k-r , 1 
k-r+1 2 



TV— >00 1 



(13) 



PoS > lim 



Pv/n-l/k 1 1 
pVH+2/k 2 


n-P-iA 1 111 1 1 

[ p+2/k IT '1 ^"p+l/kj 


-2AT-r-2AH(A + r,A) 









Notice that for large n, — > 0, because r = [y/n], and A > n. Furthermore it is AH(A + r, A) — > 0, 
because AH(A + r, A) < ^, with A > n and r = 0{y/ri]. Assuming that p is a constant, the limits 
of numerator and denominator exist, and we can use them to calculate the limit of the fraction. 
Given y > 0.5 the PoS lower bound is simplified to: 

PoS > (^1 /4 + 1 (1 - Inp)) / (3/4 + 1 (1 + Inp)) 

Numerical maximization over p G (0,1) yields PoS > 1.45, for p ~ 0.18. We also searched 
computationally for r maximizing (13), for increasing n. For r = 0.27y^ wc found the values 
appearing in the table 2(b) , that indicate PoS > 1.77. One can verify that any configuration other 
than s* and s is more expensive, by definition of distances of agents from v and v{. 



5 Approximate Strong Equilibria 

Strong equilibria do not generally exist in the Facility Location game, even for unweighted agents 
on metric networks with uniform facility costs. We illustrate this by an example, over the network 
shown in Fig. 3. Consider the metric case in Figure 3. There are three unweighted agents i = 1 ..3 
situated on distinct nodes Ut of the depicted 6-nodc cycle. For any equilibrium strategy profile s 
there exist two agents with 3 = (I + 1 ) mod 3 with Ci(s) > 1, Cj(s) > |. By agreeing to open 
a facility on vertex V|, i and j would change their costs to Ci(s') = | and Cj(s') = respectively, 
each lowering their cost by at least p. Hence, this example does not have strong equilibria. 

Existence of pure equilibria for weighted agents is also an open issue. However, we were able to 
reduce to a constant the logarithmic approximation factor a known for general network design [1, 7]. 
In fact, our result is even more general, as it concerns strong equilibria. We make use of the following 
remark. 
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Figure 3: Example of non-existence of strong equilibria. 



Remark 1 // an instance of the Facility Location game does not have strong equilibria, then there 
is at least one cycle of deviations of particular coalitions that results in a circular sequence of 
configurations {s'l-L^ with = s*^. 

Given such a sequence {s'}j^i, we denote the coalition that deviates from s' to form s'^^ by Ij. 
Such a deviation causes a cost decrease of agents in Ij and possibly a cost increase of agents in 
A \ Ij. Recall that A is the set of all agents. We define two quantities, the weighted improvement 
impr(Ij) for agents in Ij and the weighted damage dam(Ij) caused by agents in Ij respectively: 




impr(Ij) = Yl y-^l^j d^Mh) -- 
We derive an approximation factor that eliminates cycles. 

Lemma 2 Let{s^^^-f with = be a cycle of configurations in a Facility Location game instance, 
caused by consecutive deviations of coalitions. The game instance has an oc- approximate strong equi- 
librium if for all such sequences a > dammax({s')]Li); where dammaxlls^jlLi) = Toaax dam(Ij)^'^'^''. 

' ' j=1...lc— 1 

Proof. If there is no a-approximate strong equilibrium we know that there is at least one cycle 
{s'l-Li such that Vj G {1, . . . ,k- l}Vi G Ij : > a. 

Because = we have that Y^=] ) ~ ^ every agent i. Then: 

It follows that daminax({s'}jLi) > oc- Hence the lemma follows by contradiction. □ 

We derive an approximation factor as an upper bound of daminax({s'}|L^ ) for any cycle. 

Theorem 2 For every a > e there exist (x- approximate strong equilibria in the Facility Location 
game with fairly allocated facility costs, even for weighted agents and general networks. 

Proof. Wc prove that damjnax({s'}jLi ) < ^ for every cycle {s'}-^] of configurations and the result 
follows from Lemma 2. Let Ij(v) be the set of agents going to v in s', but not in s'+\ and Aj{v) be 
the set of agents going to v in both s' and s'"'"^ Note that Ij = [j^^-yl][v], therefore: 
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dammax({s']jLi; 




/c-(s'+^1^ w(ij(v)) 

daminax({s'}-Li) < max Y[ 



5,v 



^ieAj (v) 



Ci(s^) 



Hence, we need only consider what happens at the worst case node. For an agent i in Aj (v) we get 
that: 

Ci(s'+M _ + wTTTm) W(Ij(v)) 



< 1 + 



It follows that: 



dam^3x({sn^=i) < max 1 + ^ ^\ ' < lim 1 + - = e 

' ),v y W(Aj(vjj J T-^oo \ ry 



□ 



Approximate strong equilibria are also approximate pure Nash equilibria, thus: 

Corrolary 2 The Facility Location game with fairly allocated facility costs and weighted agents 
has (X- approximate pure strategy Nash equilibria for every a > e. 

6 The Strong Price of Anarchy 

We derive next an upper bound on the SPoA of a-approximate strong equilibria, for the general 
Facility Location game. 

Theorem 3 For any constant cc> e, the Price of Anarchy of oc- approximate strong equilibria for 
the Facility Location game with fairly allocated facility costs, is upper bounded tightly by 0(H(n)) 
for unweighted and by O(lnW) for weighted agents, where W is the sum of weights. 

Proof. Let s and s* be the most expensive strong equilibrium and the socially optimum configu- 
ration respectively. For any facility node v G Fj* let As* (v) and Is* (v) be respectively the subsets of 
agents connected to v in s*, and connected to v in s* but not in s. We will upper bound the SPoA 
by maxvgF;;. (HieA . (v) '-i-^^^) /Cv(s*) • Because s is an a-approximate strong equilibrium, for 
every subset of Is* (v) there is at least one agent i, that is not willing to deviate to v in coordination 
with the rest agents of the subset. Let Is* (v) = {1 , . . . , |Is* (v)|} and Ig* (v) ={!,..., i}, where i is 
not willing to deviate in coordination with I^* (v) \ {i}. Then: 

Ci s < al Wid Ui,v + . j < al Wid Ui,v + . 1 ^ 

V W IL v -h Ws v) V W IL V V 
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Figure 4: Unweighted Non-metric Lower Bound for the SPoA of a-approximate strong equihbria 
(oc > e). 



I ' V - ; 

Y_ Ci(s) < a( Y_ Wid(Ui,v) + |3v Y. 



ieis*(v) 



ieis«(v) 



i=l 



The result will follow from (14), because Is*(v) C As*(v), and agents in As*(v) \ Is*(v) play v in 
s and s*. We analyze the SPoA ratio. For Wi = 1 it is W(Ig, (v)) = |Is*(v)| = i, which yields 
SPoA = 0(aH(n)). For weighted agents set W(Is*(v)) = and use x > 1 + Inx for x G (0,1) 
on (14): 

oc > 1 \ ■ . ' < > In — /■ ' = ccfl +lnW) 

^ V W{lUv])J- ^ W(I^(v)) 

□ 



Unweighted Lower Bound We refer to Fig. 4 for a tight (non-metric) lower bounding example in 
case of unweighted agents. Facility opening costs are 1 . It is clear that all agents being connected to 
node Veq is the most expensive a-approximate equilibrium, of social cost aH(n): no agent coalition 
has incentive to deviate to Vopt, which is practically the node that all agents reside on, and results 
in a social cost of virtually 1 , for e ^ (e.g. e = n^^) . 



6.1 Unweighted Agents on Metric Networks 

We prove the following for the SPoA of the unweighted metric Facility Location game: 

Theorem 4 The Price of Anarchy of e- approximate strong equilibria in the unweighted metric 
Facility Location game with fairly allocated facility costs is upper bounded by a constant. 
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V 




Figure 5: The situation examined in the proof of theorem 4. 



The proof of theorem 4 consists of several partial results that we will synthesize. Let s denote 
any a-approximate strong equilibrium and s* the socially optimum cofiguration. We will upper 
bound the SPoA again by maxvePj* (XigA * (v) '^i(s)) / Cv(s* )^ • For any facility node v G Fs* > define 
As(v) C As*(v) to be the subset of those agents that are connected to v both in s* and s. Define 
Ig*(v) = As*(v] \ As(v] to be the subset of agents that are connected to v in s* but not in s. See 
fig. 5 for an illustration of the definitions. To simplify notation we use d(ui, v) = x|, for i G As*(v). 
At first we consider the following simple case: 

Lemma 3 Let Is* (v) he a subset of misconnected agents under ex.- approximate strong equilibrium 
profile s. Let R denote agents connected to v under s, with R C A\ Ag. (v) . // under s_r no agent of 
Is*(v) has incentive to deviate tov in coordination with Is*(v), then ^tgA ,(v)''i-(^) — (1 +C'-)Cv(s*). 

Proof. If under s_r no agent i G Ig* (v) has incentive to deviate in coordination with Ig* (v) \{i} to 
V, then for every i G Is*(v) it is Ci(s) < aCi(s*), because s is an a-approximate strong equilibrium. 
Then: 

Y_ Ci(s) = Y_ '^^(s) + Y_ '^i(s) 
ieAs*(v) ieAs(v) iels.(v) 



ieAs(v) V iels*(v] / 

< (1+a)[|3v+ Y 



The latter equals (1 + a)Cv(s*). □ 

For the rest of the analysis we treat the complementary case of that described in the previous 
lemma. Assume there exists at least one agent i G Is* (v) willing to deviate to v in coordination 
with the coalition Ij* (v) under s_r. Define a minimal disagreeing subset Ig* (v) C Is* (v) to be a 
minimal subset of misconnected agents containing an agent i that would actually deviate to v with 
Ig. (v), under s_r. Then also define Js*(v) = Is*(v) \ Is*(v). Fix i G Is*(v) to be from now on the 
agent (or one of them if there are many) that would deviate to v in coordination with Is*(v). We 
call i the unstable agent of the minimal disagreeing subset Ig* (v). By definition, the following holds 
for an unstable agent i of the minimal disagreeing subset: 
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Ci(s) > a ^xt + 1^ 



v)| + |As(v) 



(15) 



The rest of the analysis consists of bounds for agents in I^, and Js* (v) separately. In particular, 
lemmas 4 and 5 describe upper bounds for these sets respectively. 

Lemma 4 Let Is*(v) he a subset of misconnected agents under an (x-approximMe strong equilibrium 
profile s. Let Ig* (v) be a minimal disagreeing subset of Is* (v) and i G I^* (v) an unstable agent. Then: 



(16) 



Proof. By minimality of Is*(v), every agent I G Is*(v) is not willing to deviate to v under s_r 
with a coalition of size |Is*(v)| — 1. Thus for every I G Is*(v): 



ci(s) < a ^xt + 1^ 



13^ 



(v)| + |As(v)|-l 



Summing over 1°* (v) yields: 



X c,{s)<a( ) 



which is upper bounded by at most as in (16). 



□ 



Lemma 5 Let Is*(v) be a subset of misconnected agents under an OL- approximate strong equilibrium 
profile s. Let Is*(v) be a minimal disagreeing subset o/Is*(v) and Js*(v) = Is*(v) \ Is*(v). Then: 



y~ Cj(s) < a 



ieJs*(v) 



X X? + Pv(h(|As*(v)|) - H(|I°,(v)| + lAs(v)l)) 



(17) 



Proof. Without loss of generality name agents j G Js*(v) by distinct indices 1 , . . . , |Js* (v)! and 
define a series of supersets of 1°* (v), as follows: I^* (v) = 1^7^ (v) U{j}. Because s is an a-approximate 
strong equilibrium, every set I^* (v) contains an agent that is not willing to deviate to v in coordi- 
nation with Ig, (v). This agent is found either in Is*(v) — {i} or in l\*(v) \ Is*(v). We can assume 
without any loss of generality that for subset I^. (v) this agent is j; otherwise we only need to 
exchange j with some agent from I^* (v) \ {i}. One easily verifies that, by definition of a minimal 
disagreeing subset, such an exchange will not affect any of our previous results up to now. Then 
we have: 

We omit |R| and sum the inequality over j G Js*(v). The result follows. □ 

The following lemma will provide a lower bound for ^jgj *(v)^j appearing in (17). Note that 
lemma 4 provides a concrete upper bound for the cost (under strategy profile s) of agents in 
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Is*(v), as a function of their connection cost in the sociahy optimum configuration s* and the 
corresponding facihty cost (3v This is not the case with lemma 5 for agents in Js*(v); the upper 
bound (17) involves socially optimum connection cost plus pv multiplied by additional terms of 
harmonic numbers. Thus we need to determine how low the socially optimum connection cost can 
be, in order to derive an upper bounding ratio for the SPoA. 

Lemma 6 Let Is* (v) be a subset of misconnected agents under oc- approximate strong equilibrium 
profile s. Let Is*(v) be a minimal disagreeing subset o/Is*(v) and Js*(v) = Is*(v) \ Is*(v). Then: 



Z -^(^^^lM_M _ cx(h(|A,*(v)|) - H(rcxrl))), r = \lUv]\ + |As(v)| (19) 



Proof. Let i be the fixed unstable agent of I^, (v). Note that, under strategy profile s, i docs not 
have an incentive to join facility node Sj for any j G Jg* (v). Thus if j pays for Sj a share of ^ (that 
is, Sj serves Aj agents in total in s): 



Ct(s) < a fd(Ua,Sj) + y^^j < afd(ni,v) + d(uj,v) + d(uj,Sj) + ^j ^ 



Ci(s) < a (x- + X- + Cj(s)) < a I X- + X- + a I X- + — 



(20) 



.3.(V)| + |AS(V),^ y 

The latter inequality derives by usage of (18) for Cj(s), and by safely omitting |R|. Using (20) and 
the lower bound for Ci(s) given in (15), we can solve for Xy By the definition of Is*(v) in lemma 5, 
it is l4*(v)| = |I°*(v)i + j, j = 1,... ,|Js*(v)|. Then we obtain: 

Xi > max< 



{°'l+a(|lO.(v)| + |A3(v)| i + |l04v)| + |As(v)|)}' ^ ^ ' • • • ' 



Finally we sum up the latter bound over all j. Notice that x? becomes non-negative only when 
j + |I°*(v)| + |As(v)| > a(|I°*(v)| + |As(v)|). Since j + |I°*(v)| + |As(v)| is an integral value, it 
turns out that xT becomes non-negative for those values of j for which it is j + |Is*(v)| + |As(v)| > 
[adlg, (v)| + |As(v)|)]. Then, by setting r = |Is*(v)| + |As(v)|, and by summing up over all j we 
obtain the specified lower bound. □ 

Proof of Theorem 4 Now we put everything together. A lower bound on Cv(s*) is: 

ieJs*(v) iGi^*(v) ieAs(v) 
Accordingly, we obtain the following upper bound on ^^^^^ ^ Ci(s) by (16): 

Y_ Ci(s) < H T. ^^^^^ 

i6A3*(v) lGl°,{v) ieJs.(v) ieAs(v) 

iei°*(v) jeJs*(v) ieAs(v) 
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Using the latter bounds for Cv(s*) and Y.\eA *(v) '^i(^)' '^^ deduce: 

SPoA < 1 + (2a|3v + 2a ^ ^1+ Y. C5(s))/(|3v+ ^ <+ Y. 4) 



< 1 + 2a X 



|3v+ X 4 + Pv(H(|As.(v)|)-H(r)) 
jeJs*(v) 

|3v+ Y. 

iGj,.(v) 

= l+2a + 2ax — (21) 

i3v+ Y 

The latter inequality emerged firstly by substitution from (17), and secondly by removal of ^^^jo 
from numerator and denominator. We maximize (21) by taking the lower bound of the denominator, 
given in (19) and obtain: 

QD ^ 1 ^ 9 ^9 H(|As*(v)|)-H(r] 

SPoA < 1 +2a + 2a rr^ — t (22) 



Because [cxr] < (a+ l)r, and by using logarithmic bounds for the harmonic numbers: 



SPoA < 1 + 2a + 2a■ 



_J_(M}!^ _ ocln^^ + a(y +lna- 1)) 



By substituting 0.5 for y > 0.5 and a = e we can maximize numerically the resulting upper 
bounding function of y = to at most a constant. Note that, when strong equilibria exist, it 

is a = 1 . Substituting so in (22) yields an upper bounding expression similar to the one of PoS (7), 
up to constant multiplicative and additive terms. □ 



7 Open Problems 

Further investigation of the weighted game is mostly challenging: it is not known whether this 
game possesses pure equilibria. Designing a counter-example seems quite demanding, as the game 
specializes in some sense the single-sink weighted network design game studied in [7] . For this case 
pure equilibria were shown not to exist generally. Extending our analysis of the PoS and SPoA to 
the weighted metric case appears to be also non-trivial. Finally, derivation of good lower bounds 
for the PoS and the SPoA of the (non-metric) weighted case is an interesting aspect of research: 
lower bounding techniques developed in [1] do not readily apply for Facility Location. 
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